Finite size melting of spherical solid-liquid aluminium interfaces 
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We have investigated the melting of nano-sized cone shaped aluminium needles coated with amor- 
phous carbon using transmission electron microscopy. The interface between solid and liquid alu- 
minium was found to have spherical topology. For needles with fixed apex angle, the depressed 
melting temperature of this spherical interface, with radius i?, was found to scale linearly with 
the inverse radius 1/R. However, by varying the apex angle of the needles we show that the pro- 
portionality constant between the depressed melting temperature and the inverse radius changes 
significantly. This lead us to the conclusion that the depressed melting temperature is not con- 
trolled solely by the inverse radius 1/R. Instead we found a direct relation between the depressed 
melting temperature and the ratio between the solid-liquid interface area and the molten volume. 



I. INTRODUCTION 



Thermodynamic properties of a material with re- 
duced dimensionality can be markedly different from 
the bulk properties The melting temperature Tc of 
materials confined to nano-meter size can, for exam- 
ple, be significantly different from the bulk melting tem- 
perature Tm- Since the first experimental evidence^ 
melting of nano-scale metal clusters has been a case 
studyi^^i^^^^^^^^^^^, involving a number 
of experimental techniques such as transmission elec- 
tron microscopy (TEM),^^^^^ X-ray diffraction,^'^ scan- 
ning tunneling microscopy (STM)^ and calorimetric 
methods.— iiii^ii^ For a recent review on size depen- 
dent melting see Ref. fl9|. For free standing particles a 
depressed melting temperature e = {Tm — Tc)/Tm pro- 
portional to the inverse radius is typically observed. 
However, a number of new investigations have reported a 
possibility of a non-linear relationship between e and 1/R 
when the particles are sufficiently small. ^^^^^^ These con- 
tradictory results give rise to the fundamental question: 
what is really controlling the depressed melting tempera- 
ture? Several phenomenological models address the ques- 
tion of size dependent melting temperature j^-'^i^^i^^i^^i^^i^^ 
The problem is complicated by the fact that different 
crystal faces have different surface energies and by the 
occurrence of surface premelting2ii2^i2^i22> that depends 
on the surface energies. In this paper we present an ex- 
perimental sample geometry that, as we will see, avoids 
these complications. 

We use transmission electron microscopy (TEM) to in- 
vestigate the size dependent melting point of nano-sized 
cone shaped aluminium needles as a function of the apex 
angle. The interface between solid and liquid aluminium 
was found to have a spherical topology thus the problem 
of crystal facets is eliminated. For a needle with fixed 
apex angle a we find a linear relation between the de- 
pressed melting temperature e(Tc) and the inverse radius 
1/R oi the spherical solid-liquid interface similar to the 
observations for spherical particles where e(Tc) = Vt/Rin 



which case ^7 is a size-independent constant. However, by 
varying the apex angle we show that the proportionality 
constant changes significantly. This implies that the 
depressed melting is not solely controlled by the inverse 
radius. Instead, we argue that the depressed melting is 
determined by the ratio between the solid-liquid interface 
area and the molten volume. 



II. METHODS 

Samples with cone shaped aluminium needles were fab- 
ricated, from 99.99% pure aluminium, using the so-called 
ion digging method. A suspension of fine diamond pow- 
der was dispersed on the vertical edge of a 0.1mm thin Al 
semi-disk (3mm diameter). The samples were mounted 
in an ion milling machine equipped with a nitrogen cooled 
stage. Sputtering with argon ions in the disk plane results 
in formation of thin needles at areas initially protected 
by diamond particles that were eventually sputtered com- 
pletely away. With this technique cone shaped Al needles 
with apex angles of ~20-50 degrees and diameters at the 
needle point around ~20-30 nm could be fabricated, see 
Fig.lHa) and (bl). A schematic illustration of this litho- 
graphic mask technique together with some preliminary 
results on Sn needles were given in Ref. |32|. 

In order to avoid evaporation, of aluminium, during 
the in situ heating experiment the whole surface of the 
needle samples were coated with a 50-80 nm thick amor- 
phous carbon layer in a plasma coating apparatus.^ In 
addition to protecting the needles, the carbon coating 
also helps to maintain the cone shape of the liquid phase 
during the in situ experiment. All samples were exam- 
ined in a JEOL 200CX transmission electron microscope 
(TEM) using an operating voltage of 200kV and the in 
situ melting experiment was performed by making use 
of a side-entry heating holder. The temperature was 
recorded with a thermocouple (Pt-PtRh 13%) placed in 
contact with the specimen cartridge. The errors that 
could occur in the temperature measurements, due to 
for example electron irradiation, are systematically small 
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FIG. 1: (a)-(bl) Bright-field images of needles with apex angle a = 40° and 26°, respectively. The definition of the apex 
angle is shown schematically in the bottom left of (a), (bl)-(dl) Bright-field images of the aluminium needle with apex angle 
a = 26°, for three different temperatures T = 923, 926 and 927 K, respectively. In (bl) no melting is observed while in (cl)-(dl) 
a spherical melting front is visible. (b2)-(d2) Corresponding intensity profiles along the symmetry axis of the needle, (e)-(g) 
Similar needle recorded at similar temperatures but magnified six times compared to (bl)-(dl). 



and do not disturb the real physical tendencies. ^^^^^ Con- 
ventional bright-field and dark-field imaging techniques 
were used to observe the solid-liquid interface. 

III. RESULTS 

In this paper, we are going to study the melting process 
of four Al needles with apex angles a =26, 29, 36 and 42 
degrees. However, we start by inspecting the TEM mi- 
crographs of needles with a = 40° and a = 26° recorded 
at T = 923 K before the melting initiates at T = T^, see 
Fig. [U^a) and (bl). Both needles display a high degree 
of symmetry which is further supported by the symme- 
try of the thickness fringes. These TEM images there- 
fore demonstrate that the needles have an almost perfect 
conical shape. To inspect the images more carefully it is 
useful to plot the intensity profile along the central axis 
of the needle as indicated by the vertical dashed lines 
in Fig. ID^bl). Moving from the thick part of the needle 
(top of Fig. [ijbl)) towards the needle tip, the intensity 
displays a period oscillatory behavior (Fig. [Hb2)) that 



stems from the thickness that varies along the needle. 
Eventually the intensity displays a discontinuity at the 
point between the needle tip and the carbon coat as in- 
dicated by the top horizontal dashed line in Fig. [H^bl) 
and (b2). Finally, the intensity varies slowly without any 
oscillations along the carbon coat due to its amorphous 
structure, see Fig. [Hb2). 

Now we turn to discuss the finite size melting of the 
needles in the temperature range Ti < T < Tm- In 
Fig. [Hcl) and (dl) we show the same needle as in 
Fig.IUbl) but recorded at T = 926 K and T = 927 K, re- 
spectively. In both cases the tip of the needle is molten. 
This fact is also visualized by the intensity profiles in 
Fig. [TJc2) and (d2). The discontinuity between the nee- 
dle tip and the carbon coat is no longer visible and an- 
other discontinuity appears at a thicker part of the nee- 
dle, see Fig. [Hc2) and (d2). We interpret this in terms 
of a melting process where the top part of the needle is 
molten and the discontinuity of the intensity profile now 
marks the interface between solid and liquid aluminium. 
As neither the liquid aluminium nor the amorphous car- 
bon coat have any crystal structure we can not easily 
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FIG. 2: (a) Temperature dependence of the cap ratio C1Z{T) 
for the needles with apex angles a = 26°, a = 36°, and 
a = 42°. The dashed lines are fits to a constant, (b) The cap 
ratio C7l{a) as a function of the apex angle a. The dashed 
line is explained in the text. 



distinguish these two elements from our bright-field im- 
ages. This explains why the discontinuity that appears 
at the needle tip for T < Ti disappears for T > Ti. 

Next, we are going to discuss the topology of the solid- 
liquid interface. A careful inspection of the image shown 
in Fig. [T](c) reveals that the first dark fringe, in contrast 
to the other dark thickness fringes, has a circular shape. 
This fringe indicates the position of the solid-liquid inter- 
face and it strongly suggests that it has a spherical topol- 
ogy. This was further confirmed by bright-field images 
taken with six times higher magnification shown in Fig. 
[TJe)-(g). As the temperature increases from T = 926 K 
in Fig. (He) to T = 928 K in Fig. [^g), the solid- 
liquid interface area becomes larger as the melting front 
moves towards the thicker part of the needle. However, 
the solid-liquid interface remains spherical at all the ob- 
served temperatures. Finally we emphasize that during 
alternating heating and cooling cycles no hysteresis was 
observed between melting and solidification - indicating 
that the system is in a true equilibrium at all tempera- 
tures. Furthermore, the system does not seem to favor 
surface premelting since no precursor of the melting could 
be detected. 

We are now going to describe the interface between 
the solid and liquid phases in greater detail. Since the 
solid-liquid interface forms a spherical cap at each tem- 
perature Ti < T < T^, one can extract a cap height 
h{T) and a spherical radius R{T) directly from the TEM 
micrographs as shown in Fig. [TJg). In Fig. [2](a) we show 
the cap ratio C7^(T) = h{T)/R{T) as a function of the 
reduced temperature e(T) = {Tm — T)/Tm for the needles 
with a = 26°, 36°, and 42°. This reveals the following im- 
portant facts: (i) the cap ratio CIZ is essentially tempera- 
ture independent for each of the needles and (ii) the cap 
ratio CIZ decreases systematically with increasing apex 
angle a. This implies that the cap geometry remains un- 
changed as the melting front moves towards the thicker 
part of a needle with fixed a. However, the cap ratio CIZ 
varies systematically as a function of apex angle a. For 
each of the needles C7l{a) was extracted by constant fits 
to the data, see dashed lines in Fig. [2](a). The values of 
C1Z{a), given in table [H are shown in Fig. [2](b). 
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FIG. 3: (a) The depressed melting temperature e{Tc) = {Tm — 
Tc)/Tm as a function of the inverse radius for needles with 
three different apex angles a = 26° , 36° , and 42° . The insert 
compares, in double logarithmic scales, the depressed melting 
temperature e{Tc) of a needle with a = 26° and aluminium 
particles on a substrate (from Ref . Til) . (b) The melting rate 
Q{a) as a function of apex angle a. The solid line is explained 
in the text. 



We are now going to describe the size dependent melt- 
ing of the aluminium needles on a more quantitative level. 
Once the melting process is initiated for > a spher- 
ical solid-liquid interface with a radius is formed. In- 
creasing (decreasing) the temperature T' leads to a cor- 
responding increase (decrease) in the radius R' of the 
spherical interface. Thus the temperature associated 
with R' is in fact the melting temperature Tc{R) of the 
spherical solid-liquid interface. By controlling the tem- 
perature T and measuring R on the TEM micrographs, 
it is possible to track the size dependent melting tem- 
perature Tc{R) as a function of R. Figure [3] displays the 
reduced melting temperature e(Tc) = {Tm — Tc)/Tm as 
a function of the inverse radius 1/R. As for spherical 
metal particles there seems to be a linear relation be- 
tween the depressed melting temperature e(Tc) and the 
inverse radius 1/R^ i.e.. 



R{Tc) 



(1) 



where Vt{a) is a constant that depends on the apex an- 
gle a. For a comparison, the inset of Fig. [3ja) shows 
the depressed melting temperature of aluminium parti- 
cles supported by a substrate^^ and the needle shown in 
Fig. m^bl). The shared linearity might however be a co- 
incidence since the proportionality constant Vt seems to 
depend strongly on the apex angle a. By linear fits to 
the data, see dashed lines in Figl3]^a), Vt was extracted 
for all the studied needles. The result of these fits, shown 
in Fig. [HJb) and table [U reveals that Vt{a) increases sig- 
nificantly when the apex angle a is increased from 26° to 
42°. This implies that for a given radius R the reduced 
melting temperature e(Tc) may vary dramatically. Thus 
the melting of the spherical interface is not controlled by 
the local thickness of the needle. This conclusion is fur- 
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ther supported directly from the TEM images since the 
melting front does not follow the contours of the thick- 
ness fringes. 



IV. DISCUSSION 

We now turn to discuss the melting process from sim- 
ple geometric considerations assuming that the needles 
have a perfect conical shape. The cap ratio ClZ{a) = 
h{T)/R{T) relates, in this case, to the apex angle a by 



cn{a) 



KT) 
R{T) 



1 -sin(a/2 + 6>). 



(2) 



where 9 is the contact angle between the carbon coat 
and the solid-liquid interface. The fact that the cap ra- 
tio C1Z{a) is approximately temperature independent for 
a fixed apex angle a implies that remains constant as 
the melting front moves towards the thicker part of the 
needle. Assuming that we know a with a reasonable de- 
gree of precision (±1.0°) one can thus derive the contact 
angle from a constant fit to the data as shown by dashed 
lines in Fig. [2l We find that < 1.2° is systematically 
small {6/a <C 1) for all the studied needles, see table H 
In fact is smaller or comparable to the estimated error 
of a, it is therefore reasonable to assume that ^ ~ 0. The 
solid line in Fig. [2](b) shows ClZ{a) calculated from Eq. [2] 
with 6> = 0. 

Next we are going to discuss the ratio between the 
molten volume and the solid-liquid interface area. For a 
given temperature Ti < T < Tm we observed a spherical 
solid-liquid interface with radius R. Now, if the temper- 
ature increases to = T + dT, a part of the solid volume 
is going to melt. This leads to a larger solid-liquid inter- 
face with radius R' = R -\- dR. The molten volume Vm 
therefore increases to V' =Vm -\- dVm where 



dVrr. 



Jcos^ /a\ 
\"tan V2"/ 



■ Cn^{a) [3 - C7^(a)] \ nR^dR. 



(3) 

Naturally, the solid-liquid interface area Scap increases as 
weh to S'^^p = Scap + dScap with 



dScap = Cn{a)4.7rRdR 
and the liquid-coat area Sic increases with 



dSic = \ I 



tan^ 



2nRdR. 



(4) 



(5) 



We are now going to assume that the condition for melt- 
ing is that the energy gained upon melting a solid volume 
dVm has to be balanced with the increased interface area 



dS, 



cap 



dSic Hence 

K{Tc)dVm = 



IsldScap + (Tsc - llc)dSic 



(6) 



where A(Tc) ^ Lm{Tm-Tc)/Tm = Lme{Tc) is the latent 
heat per unit volume, is the latent heat of fusion. 
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FIG. 4: Reduced melting temperature e(Tc) as a function of 
surface area to molten volume ratio A{a, R) = J^{a) / R where 
jF(a) is given in the text. 



and (75c, 7k) is the solid-liquid (solid-coat, liquid- 
coat) interface energy. 

Now since the contact angle between the carbon coat 
and the solid-liquid interface is approximately zero we 
can assume, via Young's relation {jic = jsc + 7sZ cos(7r + 
6>)), that 75/ ^ 7sc - lie 

From Eqs. ([3]), (|4j), ([5|), and (|6|) we now have 



(7) 



where A{a^ R) = 
volume ratio and 



Hoc) 



LmR 

J^{a)/R is the surface area to molten 



4C7^(a) +2./1 



cos^ 
tan^ 



(f) 



(f) -cn^ia) [3-Cn{a) 



(8) 



is a geometrical factor depending only on the apex angle 
a. In table [U J^{a) is evaluated for each of the needles. 

The depressed melting temperature e{Tc) for a needle 
with fixed apex angle a is therefore, in this simple model, 
proportional to the inverse radius 1/i? as observed by 
the experiment. Furthermore the model predicts that 
the proportionality constant Q{a) = J^{a)jsi/ Lm where 
Lm ~ 867 MJ/m^ (Ref. 36). To compare the model with 
the data, the only adjustable parameter is the solid-liquid 
interface energy 7^/. The solid line in Fig. ^h) shows 
rt{a) as predicted by the model with 7^/ = 55 mJ/m^. 

The reasonable agreement between the model and the 
data indicates that the a-dependence of the proportion- 
ality constant ft{a) should be understood directly from 
the surface area to molten volume ratio A{a). To visu- 
alize this point, we show in Fig. [H the depressed melting 
temperature e{Tc) as a function of A{a^ R) for all the 
studied needles. Now all the data points collapse onto a 
single line with the slope jsi/Lm- Our results therefore 
strongly suggest that the melting process is governed by 
the surface area to molten volume ratio and not by the 
local size of the system. Notice that while this might 
seem as an obvious conclusion from a theoretical point of 
view it is not trivial to demonstrate it experimentally. 

By varying the apex angle a we succeeded to change 
the surface area to molten volume ratio A{a^ R) in a sys- 
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TABLE I: Cap ratio C7l{a), contact angle 0, proportionality 
constant ^{a), and the geometric factor J^{a) for the four 
studied needles. The cap ratio ClZ{a) and the proportionality 
constant Q{a) were derived from fits to the data shown by 
dashed lines in Fig. [2] and Fig. (Sja). Knowing ClZ{a) the 
contact angle can be calculated from Eq. (|2]) while T{a) 
can be evaluated directly from Eq. dS]). Notice that for one 
needle a = 36° the contact angle has a small negative value 
which is an unphysical result. This negative sign is related to 
the experimental uncertainty involved in determining 0. 



a 


26° 


29° 


36° 


42° 


cn{a) 


0.761(4) 


0.745(7) 


0.709(3) 


0.609(8) 


[nm] 


0.23(1) 


0.24(3) 


0.281(7) 


0.35(2) 


e/a 


0.032 


0.019 


-0.030 


0.048 




4.08 


4.34 


4.91 


5.36 



tematic fashion. This is shown in Fig. [2] - as the apex 
angle a increases the cap ratio ClZ{a) decreases, which 
in turn leads to a corresponding increase of the surface 
area to molten volume ratio A{a^ R) ~ ^(o^), see table HI 
Thus we managed to vary A{a^ R) without necessarily to 
change the radius R as opposed to the case of spherical 
particles. The advantage of the needle system is therefore 
that we could show that the depressed melting tempera- 
ture e(Tc) scales directly with the surface area to molten 
volume ratio A{a^ R) by varying both a and R. 

We now return to the discussion of the finite size melt- 
ing of metal particles on a substrate. The literature con- 
tains examples of both linear^' and non-linear— relation- 
ships between the depressed melting temperature e(Tc) 
and the inverse radius. Our results suggest that a linear 
relation is observed only if the surface area to molten vol- 
ume ratio scales with 1/R. This is of course the case for 
free standing spherical particles. However, for particles 
deposited on a substrate, that might not always be the 
case. Recently J. Murai et alS^ showed that the contact 
angle of Sn and Bi particles, supported by a substrate, 
is size dependent for < 20 nm. This indicates a more 
complicated surface area to molten volume ratio. 

We close the discussion by commenting on the ex- 
tracted value of the solid-liquid interface energy = 
55 mJ/m^ of aluminium. Previous experiments on alu- 



minium particles deposited on a substrate suggested val- 
ues of ^si that is a factor of two larger. ^^^^^^^^ These 
values should, however, be compared cautiously because 
(i) different models have been used and (ii) the shape 
of metal particles on a substrate may not be perfectly 
spherical due to the presence of the substrate i.e. the ex- 
act surface area to molten volume ratio is typically not 
specified. Finally, (in) we stress that in our model we as- 
sumed that ^si ~ Isc — lic However, if we instead assume 
Isi > Isc - lie, i-e. jsc ^ the model would capture 
the data equally well albeit 7^^ ~ 180 mJ/m^. The ex- 
act value of jsi therefore rely strongly on the relatively 
unknown interface energies of the carbon coat. 

V. CONCLUSION 

In summary, we have reported on the melting of alu- 
minium when confined to nano-sized needles. The melt- 
ing initiates at the needle tip and an equilibrium solid- 
liquid interface could be followed as the melting front 
moves towards the thicker part of the needle with increas- 
ing temperatures. The solid-liquid interface has a spher- 
ical topology i.e. the interface forms a spherical cap. For 
a needle with fixed apex angle the melting is found to be 
proportional to the inverse radius of the spherical cap. 
However, by varying the apex angle we showed that the 
proportionality constant depends on the ratio between 
the solid-liquid interface area and the molten volume. 
This lead us to conclude that the finite size melting of 
aluminium is controlled merely by the surface area to 
molten volume ratio rather than the actual local size of 
the system. 
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